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Abstract 

This work is a continuation of paper [13] where the Boltzmann 
weights for the N-state integrable spin model on the cubic lattice has 
been obtained only numerically. In this paper we present the analyti- 
cal formulae for this model in a particular case. Here the Boltzmann 
weights depend on six free parameters including the elliptic modulus. 
The obtained solution allows to construct a two-parametric family of 
the commuting two-layer transfer matrices. Presented model is ex- 
pected to be simpler for a further investigation in comparison with a 
more general model mentioned above. 
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1. Introduction 



In paper [13] a new family of N-state integrable three dimensional two-layer 
models on cubic lattice has been formulated. The weight functions of this 
model have the composite structure. Namely, they consist of eight more ele- 
mentary weights which depend on eight indices and four parameters x. These 
weights are a generalization used in [11] of the weights for Baxter- Bazhanov 
model [7]-[8]. They have the "Body-Centered- Cube" (BCC) structure firstly 
introduced by Baxter in [4] (more detailed [5]) for Zamolodchikov's model 
with N=2 ([l]-[2]). In papers [11]-[12] the idea to use a pair of so-called 
modified tetrahedron equations instead of a single tetrahedron one has been 
proposed. Paper [13] is a further generalization of this idea. Contrary to [11]- 
[12] the elementary weights should satisfy only one of two sets of modified 
tetrahedron equations. It allows to find a wider solution which provides the 
integrability of the above mentioned composite model. It may be parameter- 
ized by the elliptic functions which depend on the "angle- like" parameters 
and one more parameter k - the elliptic modulus. In these terms each elemen- 
tary weight depends on three "angle- like" parameters and modulus k instead 
of the four parameters x. Any two weights corresponding to the neighbouring 
cubes in the lattice have one coinciding "angle-like" parameter and opposite 
elliptic modulae: k and —k. Rather cumbersome analysis has shown us that 
the solution of the tetrahedron equations for the composite model has twelve 
free parameters including elliptic modulus. This result has been obtained 
only numerically. We have no analytical formulae for this general case. 

Our aim in this paper is to consider a simple particular case which could 
be more useful for a calculation of the statistical sum and other physical 
quantities. At the end of Section 6 of [13] an example of such a model has 
been considered but it seems to be trivial. In this paper we present another 
particular case which is probably more substantial. 

The paper is organized as follows: In Section 2 we remind the reader of 
the basic formulae and notations used in papers [9]-[13]. In Section 3 we 
give all necessary constraints and their solution for one composite weight W. 
Also, we mention one interesting property of this solution for our particular 
case: Z-invariance for the composite weights [6]. This property works only 
in two-layer level for the general case. In Section 4 we discuss the necessary 
conditions for the existence of the intertwining composite weights for W and 
W'. The resulting formulae are collected in the Appendix. 
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2. Formulation of the model 



We would like to recall the basic notations used in [9]- [13]. First of all let us 

remind the reader that we shall consider a spin model on the cubic lattice. In 
the Baxter-Bazhanov case the Boltzmann weights corresponding to the cells 
of this lattice are the same. This model may be called homogeneous. Each 
Boltzmann weight depends on eight indices and parameters x (see Fig.l): 
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= W{a\eJ,g\b,c,d\h;{x}) 



Fig. 1 

Up to some face factors we have 
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W{a\e,f,g\b, c,d\h; {x}) = 

w(,T3, Xi3, xi\d, h + a)w{x4, X24, X2\a, g + cr) 
w{x6,X5e,X5\e,c+ a)w{x8/uJ,X78,X7\f,b + a) Jg 



(2.1) 



where the lower index "0" after the curly brackets implies that the expres- 
sion in the curly brackets is divided by itself with all exterior spin variables 
equated to zero. The parameters Xi and Xij satisfy the Fermat condition: 
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(2.2) 



In (2.1) we used the following notations: 
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/ and k are elements of Z^. 

Below we shall use the following normalization conditions: 

X3 — X4 — Xq — Xs — 1. (2.5) 

In fact, each weight W{a\e, f, g\b,c,d\h) depends only on four continuous 
parameters Xi, X2, X5, X7 which in the case of the Baxter-Bazhanov model are 
connected with each other by the following constraint: 

Xi X2 = X5 x-j. (2.6) 

The weights functions of such form satisfy the tetrahedron equation ([9]- [10]): 



^ W{a4\c2, ci, C3I61, 63, 62!^^) W'{ci\b2, aa, 6i|c4, d, C6I64) 

d 

X W"{bi\d, C4, cslas, 63, b^c^) W"'{d\b2, 64, &3IC5, C2, Celoi) 
= ^ W'"{bi\ci,C4,C3\a2,a4:,a3\d) W" {ci\b2, as, a4\d, C2, celai) 

d 

X W'{a4\c2, d, C3\a2, 63, oilcs) W{d\ai, 03, a2|c4, C5, C6|64)- (2.7) 

where weights W, W", W" depend on their own sets of four parameters 
{x'i}, Wi}, Wi'} , accordingly. 

We write (2.7) more briefly in the shorthand notations: 

W W W" W" = W" W" W W. (2.8) 

Below we shall not demand the condition (2.6). So, each Boltzmann 
weight depends on four independent variables x. 

Authors of [11]-[12] have proposed to consider a pair of the tetrahedron- 
like equations instead of a single one (2.8): 

W W W" W" = W" W" W W (2.9) 

and 

W W W" W" = W" W" W W (2.10) 

where W depends on generally speaking some new variables Xi. Equations 
of the form (2.9) or (2.10) are called modified tetrahedron equations. This 
idea allows the commutativity of two-layer transfer matrices for the cubic 
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lattice with a "chess" structure. The composite Boltzmann weight which 
corresponds to the 2x2x2 cube of this lattice satisfies the tetrahedron 
equation. 

In our previous paper [13] we have generalized this idea in a following 
way. We have considered the composite weight W which consists of eight 
generally speaking different elementary weights of the form (2.1) as shown 
in Fig. 2: 




Fig. 2 A composite weight W. 



W^Yl W{{xa})W{{x,})Wi{x,})Wi{x,}) 

W{{x,})W{{xf})W{{x,})W{{xH}) , (2.11) 

where sum J2 is implied to be over one internal index. 

In order to satisfy the tetrahedron equation for such composite weights 
we have to consider sixteen modified tetrahedron equations for elementary 
weights of the form (ref.(6.3) of [13]): 

W W W" W" = W" W" W W. (2.12) 

We would like to note that there is no such pair of equations among them 
which would be "conjugated" to each other as a pair of (2.9) and (2.10). If it 
was so we would be forced to restrict ourselves to the model proposed in [11], 
[12] . Unfortunately, this system of sixteen modified tetrahedron equations 
appears to be too complicated to be solved analytically. We have succeeded in 
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finding the solution to this system only numerically which depends on twelve 
free parameters. Also, a particular case of tliis model has been considered at 
the end of Section 6 of [13]. In particular, the opposite elementary weights 
within the composite one are connected with each other by the inversion for 
that case. 

Now we shall consider another particular case with the following con- 
straints on the opposite elementary weights: 

Wa = Wh, = We, = Wf, Wa = Wg. (2.13) 

The same relations will be implied to be valid for another composite 
weights W, yV" yV'" . Now the system of the sixteen modified tetrahedron 
equations mentioned above reduces to the following system of the eight ones: 

WgW'gW'^W'g = w^'wlw'gWg, ^ w ^ W ^ w ' = W ^ w ^' W ^ w g , 
WeW'^W'/W'f = w'/'W"w^We, WeW'gWlwii' = W"j[wi[w'gWe, 

WhW)w'lw'l = wll'wlw'fWh, WhWlw'gW'i' = w"lw'^W'f^Wh. 

(2.14) 

So, our problem has been reduced to finding the solution to this system 
for 32 elementary weights which depend on their own sets of four parameters 
X (see (2.1)). Below we shall consider this problem in more details. 

It is more convenient to solve the system of modified equations (2.14) 
in a number of steps. First of all let us note that four modified equations 
among (2.14) have the form of (2.9) while four residual ones are of the form 
(2.10). So, our first step is to solve only one modified equation such as (2.9) 
or (2.10). Then we are going to resolve those constraints which connect with 
each other the parameters only of one composite weight W. The same can 
be done for others W", W" ,W"". The next step will be to find the necessary 
conditions on the parameters of the composite weights W, W which provide 
the existence of the intertwining weights W and W". Our last step reduces 
to applying the solution obtained in the first step to the modified equations 
(2.14) and to finding the elementary components for W" and W" . 
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3. Solution of the modified tetrahedron equa- 
tion 



Let us consider the modified tetrahedron equation (2.9). We need some 
additional notations for one elementary weight W (2.1) which depends on 
the four parameters Xi as was mentioned above. Namely , let us introduce 
the following notations: 

X2^x^, Xs^xf, Xr^x^, (3.1) 

^13 = 2^131 ^24 = X2i, X^Q = X^Q, Xyg = X^, (3.2) 

and taking into consideration conditions (2.5) 

Xx2, — Xi — 1, X^i — X2 — 1, X^Q — X^ — 1, X^s — X'j — 1. (3.3) 

Then we can introduce a set of the four parameters {m, Ji {i — 1,2,3)} 
instead of Xf 

2 X1X2 X-j X^ X56X78 , , 

m = ^r^, Ji = J2 = J3 = y y • (3.4) 

Note that m is connected with the elliptic modulus k introduced in [13] by 
the formula: 

l-k 

m = ;- 

1 + A; 

Also, it is convenient to introduce three more values 

h = m^Ji. (3.5) 

These notations are shghtly different from those which were used in [13]. 
These variables are different for the different weights. For the validity of the 
modified tetrahedron equations these variables should satisfy some definite 
relations which will be considered below. In order to make it more descriptive 
it is convenient to associate variables Jj with three cube's faces joined with a 
vertex "a" and /j with three opposite faces joined with vertex "h" as shown 
in Fig.3. 
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Fig. 3 

The connection with the "angle-hke" variables 9i introduced in [13] is as 
follows: 



J. ^ }_fcn{e,,k) 



dn{6i, k) 



m\cn{9i, k) + dn{9i, k) J ^ ^ 

In [13], all algebraic relations which provide the validity of the modified 
tetrahedron equation (2.9) where done. One can choose from them those 
relations which contain the parameters of W, W, W", W" only: 



Jl — J'2 ^ J2 — J'2} h — ^/J'-ii 

T' _ ^| I'll J' _ j'l J" _ J'" 
ll — -L/ ^3 , ^3 — ^3 , 'Jl — 'Jl 



and 



X24 x'^ 



v' v'" 

^24 ^24- 



(3.7) 



(3.8) 



Below we shall call all relations similar to (3.8) as relations of type II. They 
connect together parameters of all weights. We are interested only in a 
particular case: 

m = m' = m" = m!" . (3-9) 

The analysis of the equations (3.7-3.9) is a bit cumbersome, but rather 
straightforward and can be easily done by MATHEMATICA. It can be shown 
that we have two different solutions. One of them seems to be meaningless 
and we shall discuss only the second one. Latter has six free parameters 
including the modulus. We will not write down the manifest formulae. We 
would like to note only that the variables X'^ and X'^' of the weights W" 
and W"\ accordingly, can be expressed rationally through the variables Xi 
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and X- of W and W. The residual relations considered in [13] give us the 
variables for W, W ,W ,W . Below we show them only for W, residual 
variables can be obtained analogically: 



m 
Ji 

X5 



1/m, 



{k = -k) 

li = Jr 



3-1: 
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y {hm' - 1) 

(J3-1) ' 



^2 2/ T 

m [Ji 



X: 



11 



(3.10) 



We shall call the substitution W ^ W as T-transformation. If two weights 
W and W correspond to the cubes which join each other by one face as it 
is shown in Fig. 4, than we have 



m' — — , J2 — h, I2 — J2- 



m 



(3.11) 



The relations of this kind we will call as relations of type I. They connect 
variables only of the neighbouring weights in the lattice. 
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Fig. 4 

The equations (3.7) are relations of type I. They connect together the 
variables of two weights. Above we have only described how we can per- 
form the first stage. Let us repeat that obtaining the manifest formulae is 
straightforward. 



8 



4. "Internal" constraints on the composite 
weight 

In order to perform our following step we need some additional information. 
First of all, it is convenient for us to rewrite the equation (3.8) in a different 
form. Let us suppose that the equations (3.7) have already been satisfied. 
Then the equation of type II (3.8) is equivalent to the four relations of the 
"edge-like" parameters introduced in [13]: 

I / II , III n 

00 + ^2^ ^2 +0-2 = U 

01 - a'l + Gi[ + a[)" = 

^2 + — ol^ + ol^ = 

a3-a^ + < + <' = 0, (4.1) 
where Oj and Qq are determined by the following formulae: 

tanaj/2 — 



sn{ao, k)sn{ai, k) 
^ sn{aj, k)sn{ak, k) ' 



tan ao/2 = k sn{aQ, k)sn{ai, k)sn{aj, k)sn{ak, k), (4.2) 
and are " excesses" : 

6*1 + 6*2 + 6*3 1^ n (. ON 

ao = ^ /C, ar ^fJr - ao, (4.3) 

where /C is a complete elliptic integral of the first kind for the modulus k. 

Below we also need the transformation properties of weights (2.1) ob- 
tained in [8], [9], [12] for the rotation p on 7r/2 around the vertical axis and r 
-reflection. These two transformations are generating elements for a whole 
group of the cube's symmetry. Namely, weight (2.1) is invariant upon these 
transformations up to some face factors if variables Xi are changed by 
for p-rotation and XJ for r-reflection. Also, we need Xf' for the rotation A 
on 27r/3 around the a-h direction of the cube (Fig. 2): 

P ■ {a, e, /, g, 6, c, d, h} {g, c, a, 6, /, h, e, d} 
vP — 2Csfi. vP = 2^l2^24 vP — x±2iL^ = ■ 
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T ■ {«, e, /, g, b, c, d, h} {a, f, e, g, c, b, d, h} 

VT y V" V VT v VT V • 

^1 — ^1, ^2 — ^2, ^5 — ^7, ^7 — ^5, 

A : {a, e, /, g, 6, c, d, h} {a, g, e, /, d, 6, c, /i} 



-A ^ Xi Xrs (Xs^Xa) -jaA ^ X13 (^5-^2) 

1 Xs X24 (Xi-Xt)' ^2 X56(Xl-X7)' 




Now let us consider the eight modified tetrahedron equations (2.14). The 
analysis of the previous Section may be applied for each of them. So, we 
can write all relations of type 1 such as (3.7) and all relations of type 11 such 
as (4.1). Combining together these relations one can obtain those equations 
which are pure "internal" constraints on the variables of composite weights. 

me = nif — rUg — = m, 



Jgl — Jfl: Jg2 — Je2, Jg3 — Jh3: 

Jel — Jhl, Jf2 — Jh2, Je3 — J f3 (4.5) 

and 

go + eo + /o + /io = 

91 + /i = ei + hi 

92 + 62 = /2 + h2 

93 + h^es + U (4.6) 

where {Jei, Iei,e^}, {Jfi, Ifi, f^}, {Jgi,Ig„g^} and {JuJuuh^} are variables 
Ji,Ii defined by (3.4-3.5) and "edge-like" parameters defined by (4.2) for 
weights VFe, M^/, Wg and Wh accordingly . The analogical relations to (4.5- 
4.6) should be valid for other composite weights W", W", W". This situation 
is shown in Fig. 5. The last four relations of type 11 may be replaced by four 
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relations of variables Xe, Xf, Xg, Xh for weights We, Wf, Wg, Wh accordingly: 

Xel3 Xf2A XgYS Xh56 _ Xe2A -^/13 ^g56 XhlS _ 



Xe2A -^/13 -^556 ^ 
Xe56 ^f78 ^g24 X, 



h78 
his 



Xen ^ f24 Xg78 X 



^eI8 ^fh% ^g\Z 



h56 



Xe78 ^f56 ^gl3 ^ h24 



Xe56 Xf78 Xg24 X his 



Xgi Xei 

m'^Xgi Xei 



(4.7) 
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lf2 



./2 




Fig. 5 

The consequence of equations (4.5,4.7) are the following relations which 
appear to be very useful: 



Xfl^P Xgl^ 

^el 
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Xel 
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X. 



/7 = g 
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Xel 



W?Xg^ Xel 



X^ 



hi 



pX, 



X. 



el 



<?1 



X 



Xe2 Xgi 

A; 



-h2 



Xhb — 

Xh7 



pX, 



e2 



Xg2 Xgl 

Xg5 Xgl 



Xg2 Xgl 



q Xel X, 



9l 



m2 Xg5 Xgl 



Xfi 

Xf2 
Xf5 
Xf7 

Xhi 

Xh2 
Xh5 
Xh7 



^pX, 



= q, 



w?pX, 



gbi 



X. 



P Xel, 
Xe2 



q 



X. 



S2 



p m^Xe2 
Xel 



X. 



95 



X. 



fl2 



X. 



95 



(4.8) 



where p and q are some new parameters. One can see that substitution of 
these expressions into equations (4.5) gives us a rational constraint on p and 
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q. After one of the relations (4.7) has been satisfied we have two different 
solutions for p and consequently q. For one of these solutions all three residual 
relations (4.7) are valid automatically and we have a six-parametric solution. 
For another one, three residual relations (4.7) give us one additional quadratic 
constraint and we have five parametric solution. We should note that this 
solution contains the Baxtcr-Bazhanov [7]- [8] model and Zamolodchikov's one 
[l]-[2] for = 2 as a particular case. Unfortunately, we have not succeeded 
in achieving a final result in that case. Below we will consider only the first 
case which is much easier . In this case the following expressions for p and q 
can be obtained: 



where X^^ and X^- can be defined by applying the formulae (4.4) to We and 
Wg. It is easy to obtain a manifest formulae for variables Xfi,X fi and Xhi,Xhi 
through the six independent variables (for example, Xgi, Xg2, Xg^, X^i, Xg2 
and m) using (3.10), (4.4), (4.8) and (4.9). 

To conclude this Section we would like to mention a remarkable fact 
concerning our solution. Namely, let us consider two neighbouring weights 
within the composite weight W (see Fig.l),for example. We and Wd which is 
now equal to Wg . We can find the intertwining weights W" and W'" for this 
pair. In order to do so we should substitute Xei and Xgi into the equations 
(4. 4), (3. 7) instead of Xi and X- accordingly and resolve them with respect 
to X'/ and X^", as it was described in the previous Section. We can follow 
the same procedure for pairs Wa = Wh and Wf , We = Wf and Wh , Wg and 
Wb = We and obtain the F/' and Y/" , and , 1//' and V^" accordingly. 
It is interesting to note the following relations : 



One can check that the similar situation takes place for all pairs of the neigh- 
bouring elementary weights inside the composite weight. This property is 
nothing else but Z-invariance generalized on the three dimensional case ([6]). 
One should note that this property works for a general model (beyond the 
condition (2.13)) only in the two-layer level and breaks within the composite 
weight. We did not demand the validity of Z-invariance for the model (2.13). 
That is why we were very surprised to observe it. 



P = 



{X^l Xgg) Xg2 

i^eb ~ ^92) Xe2XgiXg5 ' 



Q = 



(X'e^-X^,) rn?Xe^Xg, 



(4.9) 




(4.10) 
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5. Existence of the intertwining weights 

Let us consider the question about the existence of the intertwining com- 
posite weights W" and W'" for two original weights W and W' which can 
be associated with the two cubes joining each other by the face as shown in 
Fig.6: 




Fig. 6 

As was mentioned above each of the weights W and W" has six free 
parameters (including modulus). Besides, we have three relations of type I: 

m = m', Jg2 = Jg2: Jf2 = J'f2- (5.1) 

It seems to be natural that there are no another constraints on the parameters 
of W and W. If it was so we would have three parametric family of the 
commuting two-layer transfer matrices. But accurate analysis has shown us 
that it is not the case. There is one more constraint which is necessary for the 
existence of W" and W'" . So, the commuting family is only two- parametric. 

After the analysis of the previous Section has been performed it is enough 
to resolve the following set of twelve equations of type I: 



Jgl 


- T'" 


Jg2 - 


J' 


JgZ - 






4i = 














Jel 


- T'" 


Je3 = 






- 1 / T'" 

— ^1 ^e3i 




J'eS 




J/2 = 


Jl 


J" - 


T" 


(5.2) 



and four relations of type II: 

^g2A ^g2A — ^g2A ^g2A-i 
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Xe24: 


Y" 


— Y' 

— ^e24 


Y'" 


Xf24 


Y" 
^h2A 


— Y' 

— ^h24 


Y'" 


Xh24 


Y" 


= x:'^2A 


Y'" 
^h2A- 



(5.3) 

Below we describe our final result. Our choice of independent variables is as 
follows: 

modulus m, five parameters from W: Xgx^Xgi^ Xg^, Xei,Xe2 and two param- 
eters from W': X'g^,Xg2- In order not to encumber the text we have collected 
the resulting formulae for X'^^ in the Appendix. 

Variables X^^ can be obtained from the formulae (3.10) and definitions 
(3.4) applied to weight W^. Expressions for Xj,j,Xj,j and X^j,X;jj can be 
extracted by substituting already known values into the formulae (4.8) for 
the composite weight W. Also, it is necessary to use the following expressions 
for p' and q': 

p'=f .' = ff^. (5-4) 

Now we know all variables for weights W and W" and use the analysis of 
Section 3 for finding all intertwining elementary weights from (2.14) using a 
known solution of the equations (3.7-3.8). It is remarkable to note that all 
"internal" constraints described in Section 4 are satisfied automatically. Let 
us remind the reader that we have solved all necessary equations for A^'th 
powers of the original variables x. Now we need to choose the right powers 
of u when extracting the A^'th roots. We have done it for a general situation 
and a particular case considered in Section 6 of [13]. This experience tells us 
that it can be done. Moreover, as a rule there is some arbitrariness in this 
procedure. 



6. Conclusion 

In this paper we have considered a particular case of the general two-layer 
integrable model proposed in [13]. In fact, we have discussed only one out 
of two possible solutions. This solution has five free parameters and one 
modulus which is the same for all composite weights. It is interesting to note 
that for this solution the so-called Z-invariance has been restored inside each 
composite weight while for the general model this property takes place only 
on the two-layer level. Finding of the intertwiners for a pair of the composite 
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weights W and W" appeared to be unexpectedly difficult. In spite of our 
expectation the commuting family of the transfer matrices has only two free 
parameters (not three). Unfortunately, we have not clarified completely the 
meaning of an additional constraint. The model discussed above is simple 
enough and at the same time it has enough free parameters. So, we hope 
that this model may be useful for study of it's thermodynamic properties. 
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Here we present the resulting formulae for X'^^. These formulae being written 
in terms of the independent variables are rather cumbersome. To simplify 
them we need more notations: 



7. Appendix 



Xg5 Xgi3 Xg24 {1- Igs): 



Xg5 Xgl?, Xg24 (1 — Jgs), 
Xg2 Aei3 Xe24 (l-/e3), 




(7.1) 



Note that 



Xgl _ Vg X. 



el 



Ve ^sl _ y'g 



(7.2) 



Xgl Zg Xf 



el 



Also, let us introduce : 



a 



TV? Xg5 yg - Zg 

Xg5 Ze — Xel He 



(7.3) 
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and 



wi = X'g^ Zg - Xg^ yg, 

W2 — Xgi Xgi Xg2 — Xe2 ^g5j 

W3 = Xg5 y'g - Xgi Zg, 

W4 = Xel Xg2 y'g — Xg5 Xe2 z'g, 

W5 = ^e24 {Xg2 - Xg^) {Xg2 - Xe2 Xg^) W2 + 

+ {X'g2 — Xe2) Xg5 (Xgi Xg2 Ze — m Xe2Xg^ye). (' 

After that our result looks as follows: 

-\ri C C V Q Q 

^el ^ '-'1 '-'2) ^e2 — ^g2 '-'1 '-'3, 

A(,2 rn Ag5 



where 



and 



Oi — — , 02 — — , 03 — 

ri r2 rs 



h = X'yi {X'g2 - Xg2) Vea- X'g2 wi 13, 
ri = Xg5 {X'g2 - Xg2) Ze a- Xg2 wi P, 

I2 = Wi W2 y'g - Xe2 Xg5 a, 

r2 = Xg2 {Xei - X'yi) W3 yg + wi W4, 

Is = Xe2 Xg2 Uh, (3 - {X'g2 - Xg2) W2 z^ y'g, 

rs = (m^-l) Xg2Xg,iXei-X'gi) {X'g2-Xg2)w5 + 

+ Xg2iX'g2-Xe2)ri + h. (' 
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